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ABSTRACT 

The fate of the cooling gas in the central regions of rich clusters of galaxies is not well 
understood. In one plausible scenario clouds of atomic or molecular gas are formed. 
However the mass of the cold gas, inferred from measurements of low energy X-ray ab- 
sorption, is hardly consistent with the absence of powerful CO or 21 cm lines emission 
from the cooling flow region. Among the factors which may affect the detectability of 
the cold clouds are their optical depth, shape and covering fraction. Thus, alterna- 
tive methods to determine the mass in cold clouds, which are less sensitive to these 
parameters are important. 

For the inner region of the cooling flow (e.g. within the radius of ~50-100 kpc) 
the Thomson optical depth of the hot gas in a massive cooling flow can be as large as 
~ 0.01. Assuming that the cooling time in the inner region is few times shorter than 
the life time of the cluster, the Thomson depth of the accumulated cold gas can be 
higher accordingly (if most of the gas remains in the form of clouds) . The illumination 
of the cold clouds by the X-ray emission of the hot gas should lead to the appearance 
of a 6.4 keV iron fluorescent line, with an equivalent width proportional to tt- The 
equivalent width only weakly depends on the detailed properties of the clouds, e.g. 
on the column density of individual clouds, as long as the column density is less than 
few 10 23 cm -2 . Another effect also associated exclusively with the cold gas is a flux in 
the Compton shoulder of bright X-ray emission lines. It also scales linearly with the 
Thomson optical depth of the cold gas. With the new generation of X-ray telescopes, 
combining large effective area and high spectral resolution, the mass of the cold gas 
in cooling flows (and it's distribution) can be measured. 

Key words: galaxies: clusters: general - cooling flow - intergalactic medium - X-rays: 
galaxies. 



1 INTRODUCTION 



The radiative cooling time of the gas in the central parts 



the Hubble time (|Lea et al. 1973; silk 1976 


; Cowie and Bin- 


ney 1977 




Fabian and Nulsen 1977) and the gas may cool 



down below X-ray temperatures, forming a cooling flow 
(see Fabian, 1994 for review). Strong peaks in the surface 
brightness distribution, observed in many clusters of galax- 
ies (e.g. Branduardi-Raymont et al. 1981, Fabian et al. 1981, 
Canizares, Stewart & Fabian A.C. 1983, Fabricant & Goren- 
stein 1983, Jones & Forman 1984, Stewart et al. 1984), are 
usually considered as indicators of a cooling flow. Evidence 
for cool gas was also found in the spectroscopic observations, 
which revealed the presence of emission lines, characteristic 
of gas with the T e < 5 x 10 6 K gas (e.g. Canizares et al., 



1979). Detailed studies of the surface brightness distribution 
in the cooling flows lead to the conclusion, that a fraction 
of the gas drops out from the flow at different radii, such 
that the mass deposition rate is proportional to the radius: 
M{< R) xR (e.g. Thomas, Fabian and Nulsen, 1987). The 
fate of the gas which has dropped out from the flow is not 
well understood. Only small fraction of the cooling gas can 
form stars with a normal IMF. The optical and UV obser- 
vations of the central galaxy restrict the total rate of mas- 
sive star formation rate typically by ~ 5-30 M^jyear (e.g. 
Fabian et al. 1984, McNamara & O'Connel 1993, Cardiel, 
Gorgas, Aragon-Salamanca 1995, Smith et al. 1997). De- 
tection of excess absorption in several cooling flow clusters 
(e.g. White et al., 1991) seems to support the hypothesis 
that the cooling matter forms clouds of atomic or molecular 
gas. To provide effective absorption, cold clouds must have 
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a covering fraction close to unity (White et al. 1991). Recent 
analysis by Allen and Fabian (1997) of a sample of nearby 
cooling flow clusters observed with ROSAT confirmed pres- 
ence of excess absorption, although with covering fraction 
somewhat lower than unity. The de rived masses of X- ray- 
absorbing material are ~ 10 12 M flWhite et al. 199l| ) and 
are hardly consistent with the absence of powerful CO and 
HI emission (e.g. O'Dea et al. 1994, Voit and Donahue 1995). 
Thus, independent determinations of the cold gas mass are 
important in understanding the fate of the cooling gas. Geo- 
metrically small, optically thick clouds may easily escape de- 
tection. Therefore observations at the wavelengths at which 
the clouds are transparent are the best indicators of total 
mass of cold gas. From this point of view, X-rays with ener- 
gies above ~6 keV could be useful, since clouds only become 
opaque at these energies when hydrogen column density ex- 
ceeds few 10 23 -10 24 cm~ 2 . 

The Thomson depth of the cooling flow regions, calcu- 
lated for the hot X-ray emitting gas, can be of the order 
of 0.01. Assuming a cooling time in the inner region which 
is few times shorter than the Hubble time, the averaged 
Thomson depth of the gas deposited during the life time of 
the cluster can be even higher. We argue below that such 
an amount of cold gas is sufficient to produce the emission 
from the 6.4 keV fluorescent iron line with a flux well above 
the sensitivity limit of future X-ray missions. Possible pres- 
ence of the 6.4 keV line in clusters was also noted earlier 
by Vainshtein & Sunyaev (1980) and White et al. (1991). 
Another effect, associated with cold gas, which also might 
be observed for such values of Thomson optical depth is 
the Compton shoulder of the brightest X-ray emission lines. 
Note that (unlike CO and 21 cm lines or low energy absorp- 
tion) the fluxes in the 6.4 keV line and Compton shoulder 
do not depend strongly on the shape or column density Nh 
of the individual cold clouds as long as Nh is lower than 



Table 1. Equivalent width of the 6.4 keV iron fluorescent line, cal- 
culated assuming an emission spectrum of optically thin plasma 
at different temperatures, and normalized to unit Thomson opti- 
cal depth and solar abundance of iron (see equation ^) 



Temperature 
(keV) 


EW (keV) 

f°° (E K /Ef 8 I(E)dE 

71 K 

U -'^ J(6.4) 


1 


0.24 


2 


0.54 


3 


0.71 


1 


0.83 


10 


1.09 



(gj) in a more convenient form: 

A Fe \ fZ(EK/Ef*I(E)dE 



EW = 0.74 



( Are \ 
U.68 10- 5 J 



1(6.4) 



tt keV (2) 



For hard power law spectra, the equivalent width of the 
6.4 keV line can be of the order of tt keV. In the cooling 
flow region, X-ray emission of the gas is relatively soft and 
equivalent width may be lower (see however next section). 
We calculated the expected equivalent width assuming, for 
I{E), an emission spectrum of an optically thin plasma at 
different temperatures (Table 1)Q Table 1 shows that al- 
though "hotter" spectra have a clear advantage in terms of 
producing a larger equivalent width of the iron line, this 
temperature dependence is not too steep and a change of 
temperature from 10 to 1 keV results in a decrease by only 
a factor of ~ 4 in the equivalent width. 



2 EQUIVALENT WIDTH OF THE 6.4 KEV 
LINE 

Let us consider the simplest case of a uniform sphere of neu- 
tral (or molecular) gas with an X-ray continuum source at 
the center. We assume that the optical depth of the sphere 
is less than 1 for X-rays above ~ 6 keV. Absorption of 
the X-ray photons above the iron K-shell absorption edge 
(Ek = 7.1 keV) will cause (with probability Y ~ 0.3, e.g. 
Bambinek et al., 1972) emission of 6.4 keV photons. Thus, 
the equivalent width of the 6.4 keV line for the whole region 
can be calculated as (e.g. Vainshtein & Sunyaev, 1980): 



EW 



YN H A Fe f£ x a Fe (E)I(E)dE 
f(6A) 



(1) 



where Nh is the hydrogen column density of the region, Ap e 
is the iron abundance relative to hydrogen (we adopt A Fe = 
4.68 10~ 5 , Anders and Grevesse, 1989), a Fe (E) = 3.5 x 
10 -20 (Ek /E) 2 " 8 cm 2 is the photoabsorption cross section 
(from the K shell) of an iron atom, 1(E) is the incident X-ray 
spectrum (in phot / s / sm 2 / keV) . One can rewrite equation 



3 THOMSON DEPTH OF THE COOLING 
FLOW REGION 

To estimate the expected equivalent width one should mul- 
tiply the value, given in Table 1, by the Thomson depth of 
the region and the iron abundance relative to the solar value. 
Assuming that all gas, deposited in a cooling flow during a 
Hubble time tn ~ 2 x 10 10 years, forms clouds of neutral 
gas, one can write: 



T T = 1.7 X 10" 



— Y 

50kpc J 



M 



100M© /year 



(3) 



The above estimate assumes that the cold gas is uniformly 
distributed over the volume and the source of continuum 
radiation lies at the center. Multiplying eq.(^) by the factor 
from Table 1 (e.g. for T e = 3 keV see also discussion below) 
one can write for a final estimate of the equivalent width: 



EW ~ 12 



— Y 

50kpc J 



M 



IOOMq /year 



)( 



A, 



4.68 10- 5 



eV (4) 



* Here and in the subsenuent sections w 




used MEKA model 


(Mewe, Gronensr 


hild van den Oord 1985 




Mewe, Lemen & van 


den Oord 1985 




Kaastra 1992 


) as implemented in XSPEC vlO 



thin plasma. 
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Note that although the hot gas in clusters has (on average) 
an iron abundance less than solar (~ 0.3-0.5), a central con- 
centrations of heavy elements was observed in a number of 
clusters with cooling flows (e.g. Fukazawa et al., f994). Note 
also that the above estimate is insensitive to the shape of 
the clouds. As long as each individual cloud has an opti- 
cal depth (for Thomson scattering and for photoabsorption 
above 6 keV) much less than 1, a clumpy distribution (with 
an arbitrary covering fraction) will result in the same equiva- 
lent width as a uniform distribution. The radial distribution 
of the scattering media is more important. 



4 DISCUSSION 
4.1 6.4 keV line 

Let us consider the case of an isothermal cluster, i.e. whose 
emission spectrum at any place in the cluster is character- 
ized by the same value of temperature T e . Assuming that 
the mass deposition rate is proportional to the radius (e.g. 
Thomas, Fabian and Nulsen, 1987), one can see that the 
most favorable conditions for producing 6.4 keV fluorescent 
iron emission is in the very center of the cooling flow. Since 
the surface brightness of cooling flow clusters is very strongly 
peaked, then the emission from this region may dominate 
the spectrum even if the spatial resolution of the instru- 
ment does not allow one to fully spatially resolve the cool- 
ing flow region. The equivalent width of the 6.4 keV line 
reflects the total mass of iron in the region and is insensitive 
to the shape or size of the clumps, as long as each individual 
clump is optically thin for X-rays above 6 keV. On the other 
hand, the same clouds can be optically thick to low energy 
X-rays and the efficiency of low energy absorption can be 
low (recalculated per hydrogen atom). Thus, if clumps are 
dense and have a small covering fraction^], they may not 
contribute significantly to the low energy absorption, but 
still effectively produce 6.4 keV photons. Note also, that the 
resulting equivalent width is almost insensitive to the tem- 
perature of the clumps, as long as the temperature is low 
enough. 

Allen et al. (1996) identified three distant clusters: 
Zwicky 3146 (z = 0.291), Abell 1835 (z = 0.252) and 
E1455+223 as having the largest known mass deposition 
rates of 1400, 2300 and 1500 M Q /year. The equivalent 
width of the line from these clusters can be of the order of 
tens of eV. For nearby clusters, A478 (z — 0.09), has a large 
mass deposition rate of order 1000 Mq /year at R ~ 200 kpc 
(Allen et al., 1993) and the equivalent width within the 100 
kpc central region may be of the order of 10 eV. One can 
make a "more accurate" estimate of the equivalent width 
for this cluster (Fig.Q), using the results of the deprojec- 
tion analysis by Allen et al., 1993 and White et al., 1994. 
Shown in the bottom panel of Fig.jj] is the expected equiv- 
alent width of the 6.4 keV line, calculated as a function of 
the projected radius from the cluster center. In these calcu- 
lations we assumed that the gas in the cluster is isothermal 

T Discussion of the constraints on the cloud parameters, which 
can be derived considering in detail conditions in the cooling rc- 
gion(see e.g. Daines, Fabian and Thomas 1994) are beyond the 
scope of this article. 
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Figure 1. Upper panel: The assumed radial distributions of 
hot gas density n e (r) in cm -3 (thin solid line), mass deposition 
rate in Mq per year (thick solid line) and the density of cold gas 
n C old(i') in cm -3 (dashed line) accumulated for 2 X 10 10 years. 
Lower panel: The equivalent width of the 6.4 keV line as a 
function of the projected distance R from the center. For this 
figure we adopted the results on n e (r) and M of White et al. 
1994. 



with a temperature of T e ~ 3 keV (Allen et al., 1993). Of 
course, given the lack of knowledge of what is the real fate 
and distribution of the cooling gas, these calculations may 
not provide more accurate estimates of the line flux than the 
simple estimates in the previous section for the region as a 
whole (assuming a uniform distribution of the cold gas) . 

The total flux in the fluorescent line is defined by the 
quantity n coH (r) x J^ k (Ek / 'E) 2 ' s £(E,r)dE (where n cold (r) 
is the density of cold gas and £(E, r) is the radiation density) 
integrated over the volume occupied by the cold gas. Note 
that only that part of the radiation density spectrum above 
6 keV determines the flux and equivalent width of the line. 
Although the emission measure weighted temperature is low 
in the cooling flow region, compared to the outer hot regions, 
it is likely that cooling flows are inhomogeneous and hot gas 
is present even at small radii (e.g. Fabian and Nulsen, 1977, 
Nulsen 1986, see Fabian, 1994 for review) and above 6 keV 
the contribution to the radiation density from this hot phase 
may be dominant (see e.g. Allen, Fabian & Kneib 1996). One 
can estimate the "effective" temperature, i.e. the value of T e 
which should be used for calculating the line flux (see Table 
1), using the high energy (above 6 keV) part of the observed 
spectrum towards the center of the cluster, since the shape 
of the radiation density spectrum at the cluster center and 
the observed spectrum along this line of sight coincide for a 
symmetric cluster. Note also, that if an AGN (having hard 
X-ray spectrum) is present in the nucleus of the central 
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galaxy, then the equivalent width of the fluorescent line will 
be appropriately higher. 

Assuming that the trend Mq oc R holds down to small 
radii and cold clouds remain in the region, where they were 
deposited, then equivalent width will be highest in the very 
central region of the cooling flow. As noted above if the ra- 
diation density (above 6 keV) is strongly peaked at the cen- 
ter of the cooling flow then the emission measure weighted 
equivalent width from the whole region (within given radius) 
can correspondingly be higher than the equivalent width cal- 
culated for given projection radius. 

With present instruments (e.g. ASCA) the detection 
of the 6.4 keV line from the cooling flow regions is ques- 
tionable. The expected flux in the 6.4 keV line is relatively 
weak compared to that of the 6.7 keV line complex and the 
combination of energy (and angular) resolution and sensi- 
tivity may not be sufficient to detect EW ~ 10 eV fea- 
tures from clusters. The better energy resolution of ASTRO- 
E (with bolometers), Spectrum-X-Gamma with the Bragg 
spectrometer - OXS (see e.g. Christen sen et al. 1990) and 



especia lly CO NSTELLATION /HTXS (White, Tananbaum 



& Kah i 1997) will allow one to detect this line from the 
clusters with massive cooling flows (if most of the cooling 
gas indeed forms molecular clouds). For Perseus (assuming 
an equivalent width of ~ 10 eV) about 100 line photons 
could be detected in a 100 ksec observation with the Bragg 
spectrometer of Spectrum-X-Gamma. Shown in Fig.[] is the 
expected number of photons which will be detected from 
the central 2' region of the Perseus cluster. An effective area 
of ~100 cm 2 (taking into account the peak reflectivity of 
~21%) and an energy resolution of ~5 eV were assumed. 
The intrinsic structure of the 6.4 keV line (which consists of 
two components separated by ~13 eV) was used to generate 
the figure. A few 100 ksec exposures (to measure line and 
continuum fluxes) should be sufficient to detect the 6.4 keV 
line with an equivalent width of a few eV. For the AXAF 
High Energy Transmission Grating extended nature of the 
source will significantly complicate the detection of the spec- 
tral features. CONSTELLATION with anticipated effective 
area of ~ 6000 cm -2 and an energy resolution of few eV 
will be able to detect very weak lines during short exposure 
thus allowing studies of clusters with much weaker cooling 
flows. With CONSTELLATION it should be even possible 
to detect 6.4 keV line associated with column densities of the 



gas 



10 cm der ived from low energy absorption 



cold 



studies (White et al. 1991) 



This part of the spectrum around the 6.4 keV line should 



be rela 



lvely free from contamination by bright 



emission 

lines from the hot gas. At temperatures lower than ~ 1 keV, 
excitation of the inner shells of FeXVII-FeXX ions by elec- 
tron impact may produce lines at 6.42-6.47 keV, but their 
equivalent widths should not be high, if the higher tempera- 
ture components (above T e ~ 1 keV) dominate the spectrum 
at energies of ~ 6 keV. 

Along with the 6.4 keV (K a ) line the Kp line at 7.06 
keV should be present with intensity about 10 times lower 
(j^q) than that of the K a line. The absorption edge at an 
energy of 7.1 keV should also be present. However it may be 
considerably more difficult to detect relatively wide absorp- 
tion feature compared to emission line. 



Perseus cluster, 2' circle 




6.4 6.5 
Energy, keV 



Figure 2. The number of photons, which could be detected by 
OXS of Spectrum-X-Gamma at a given energy during a 10 5 s 
exposure towards the Perseus cluster (central 2 arcmin circle). 
The value of 10 eV was adopted for the total equivalent width of 
the 6.4 keV line (which consists of two components). The energy 
resolution of the OXS is ~ 5 eV. 



4.2 Compton shoulder 

Another effect, which also scales linearly with the Thom- 
son depth of the cold gas is the Compton shoulder, associ- 
ated with the brightest emission lines. Assuming a Thomson 
depth for the cold gas of tt ~ 0.1 one can conclude that ~ 
10% of all photons (in particular, those in emission lines) 
will be scattered, forming the so-called Compton shoulder 
at the left (low energy) side of the emission lines (Fig. ||). 
Hot gas does not contribute to the formation of the Compton 
shoulder, since large thermal velocity of the electrons causes 
complete smearing of the features. On the other hand scat- 
tering of the emission lines by electrons bound in atoms and 
molecules will produce the Compton shoulder wit h distor- 
tions, specific to th e given type of atom or molecule ( Sunyaev 
Churazov 199(| ). This effect is also weakly dependent on 



the cloud parameters and in principle provides the unique 
possibility of distinguishing between the contributions of free 
(~ cold) electrons, atomic hydrogen, atomic helium or singly 
ionized helium, although presence of blends of emission lines 
(like the 6.7 keV complex) may significantly complicates the 
measurements. In principle, combined measurements of the 
iron K a line and the Compton shoulder allow one to deter- 
mine the iron abundance (tet) in the cooling flow. 



5 CONCLUSIONS 

Future X-ray missions, combining high energy resolution 
and large effective area, will be able to estimate the amount 
of cold gas deposited in massive cooling flows, by measur- 
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Figure 3. The dotted line shows the emission spectrum from an optically thin plasma (MEKA model as produced by XSPEC) at T e = 2 
keV (left) and T e = 4 keV (right), convolved with a FWHM=5 eV Gaussian. Solid line shows the effect on the spectrum of scattering 
by neutral hydrogen (for a Thomson optical depth of 0.1). The scattered ~ 6.6—6.7 keV photons have slightly lower energies, causing an 
increase in the continuum level at 6.5—6.6 keV by a factor ~ 1.5 X (Sj-))- Note that in the presence of lower temperature components 
this 6.5—6.6 keV region may also be contaminated by weaker emission lines. 



ing the equivalent width of the 6.4 keV iron fluorescent line 
or the intensity of the Compton shoulder, associated with 
brightest X-ray emission lines. Assuming that the bulk of 
the cooling gas forms cold molecular clouds, the expected 
fluxes are high enough to ensure their detection. These mea- 
surements could provide robust estimate of the gas mass, 
insensitive to the detailed shape, density or temperature of 
the cold clouds, as long as column densities of individual 
clouds are less than ~ 10 24 cm -2 . In particular, the inten- 
sity of the 6.4 keV line will be less dependent on the covering 
fraction of the clouds than either the magnitude of the low 
energy absorption or the intensity of the 21 cm or CO lines, 
both of which are more sensitive to the column densities of 
individual clouds. 

This work was supported in part by the RBRF grant 
96-02-18588. C. Jones and W. Forman acknowledge support 
from the Smithsonian Astrophysical Observatory and AXAF 
Science Center. 
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